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Sums of matrix elements of spin-dependent two-body momentum-independent interactions and 
sums of their products are calculated analytically in the basis of many-body states with given 
total spin — the states built from spin and spatial wavefunctions belonging to multidimensional 
irreducible representations of the symmetric group, unless the total spin has the maximal allowed 
value. As in the first part of the series [V. A. Yurovsky, Phys. Rev. A 91, 053601 (2015)], the 
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the Hamiltonians of non-interacting particles. The sum rules are applied to perturbative analysis of 
energy spectra and to calculation of two-body spin-dependent local correlations. 
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INTRODUCTION 

The present paper continues analyses [1] (see this work 
for a more comprehensive introduction and definitions) of 
the sum rules for many-body systems of indistinguishable 
spin-i particles. The particles can be composite, e.g., 
atoms or molecules, and the spin can be either a real 
angular momentum of the particle or a formal spin, whose 
projections are attributed to the particle’s internal states 
(e.g. hyperfine states of atoms). In the latter case, the 
spin i means that only two internal states are present in 
the system. This formal spin is not related to the real, 
physical, spin of the particles, which can be either bosons 
or fermions. 

The many-body wavefunctions are represented here, 
as well as in [1], as a sum of products of the collective 
spin and spatial functions. These functions depend on 
spin projections and coordinates, respectively, of all par¬ 
ticles and belong to multidimensional, non-Abelian rep¬ 
resentations of the symmetric group (see [2-5]), unless 
the total spin has the maximal allowed value. For spin- 
i particles, the representation is unambiguously deter¬ 
mined by the total many-body spin. This approach is 
applicable to the Hamiltonians which are separable to 
spin-independent and coordinate-independent parts. It 
differs from the conventional approach (see [6] and [1]), 
where each particle is characterized by its spin projec¬ 
tion and coordinate, and the total wavefunction is sym¬ 
metrized for bosons or antisymmetrized for fermions over 
permutations of all particles. The total many-body spin 
is undefined in this case. As follows from the Heitler’s re¬ 
sults [7] (see also [I]), an exact wavefunction for particles 
with spin-independent interactions can not be obtained 
in the approach with defined individual spin projections. 
Many-body states with defined total spin [8-15], includ¬ 
ing the collective spin and spatial wavefunctions [16-19], 
were applied to spinor quantum gases (see [20-30]). Such 
states were also proposed for implementation of permu¬ 
tation quantum computers [31]. Other kings of entangled 
states with non-trivial symmetry have been analyzed for 
quantum-degenerate gases of spin-1 [32] and spin-2 [33] 


bosons. A more comprehensive review is presented in 
Ref. [1]. 

The well-known mean-field approach (see [6, 24]), 
where the interactions between particles are replaced by 
the self-consistent field, is generally applied to states 
with defined individual spin projections and undefined 
total spins. It can be applied to certain states with de¬ 
fined total spins if the interactions are spin-independent. 
For bosons, the many-body mean-field wavefunction con¬ 
tains a single spatial orbital, determined by the Gross- 
Pitaevskii equation (see [24]). For fermions, the Hartree- 
Fock wavefunction (see [6]) contains double-occupied spa¬ 
tial orbitals. Such wavefunctions have the total spin 
S = N/2 or S' = 0 for iV spin-i bosons or fermions, re¬ 
spectively, and describe quantum-degenerate gases with 
multiple occupations of the spatial orbitals. The present 
approach describes states with arbitrary total spins and 
is suitable for non-degenerate gases. 

For spin-independent two-body interactions between 
particles, sums of matrix elements directly follow from 
the Heitler’s results [7], while the sum of the matrix 
element squared moduli was calculated in [1]. Matrix 
elements of spin-dependent interactions depend on the 
total spin projections of the coupled many-body states. 
For spin-dependent one-body interactions with external 
fields, this dependence was factorized using the Wigner- 
Eckart theorem and sums of the matrix elements and 
their squared moduli were calculated [1]. 

The present work is devoted to spin-dependent two- 
body interactions between particles. As well as in Ref. 
[1], the matrix elements are calculated in the basis of 
non-interacting particles with single occupation of spa¬ 
tial modes. The interactions are expressed in terms of ir¬ 
reducible spherical tensors and their matrix elements are 
related to ones for the maximal allowed projections of the 
total spins using the Wigner-Eckart theorem in Sec. I. 
Sums of these matrix elements and their squared moduli 
are calculated in Sec. H for zero-range spin-dependent 
interactions. The sum rules are applied to perturbative 
analysis of energy spectra in Sec. HI. Average two-body 
spin-dependent correlations are analyzed in this section 
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too. These correlations allow to distinguish between the 
many-body states with defined total spins and individ¬ 
ual spin projections. Appendices contains calculation of 
sums, used in Sec. II. 

References to equations in the previous paper [I] are 
started here from I. The present paper, as well as [1], uses 
the following notation for the universal factors. In the 

ratios of the 3j-Wigner symbols giving the de¬ 

pendence on the total spin projections, S and S' are the 
maximal and minimal total spins of the coupled states, 
q is the rank of the spherical tensor, and Sz and Sz + k 
are the spin projections for S and S', respectively. The 
factors and B], respectively, appear 

in the sum rules for the matrix elements of A and prod¬ 
ucts of the matrix elements of A and B for the maximal 
allowed spin projections. Here A and B are arbitrary op¬ 
erators, S is the maximal total spin of the coupled states 
and m is the number of changed spatial quantum num¬ 
bers (it is omitted if the factor is independent of m). If 
the sum rule contains several factors Y, they are specified 
by the subscript i. 

I. THE SPIN-PROJECTION DEPENDENCE 

Let I t) and | j,) denote two spin states of the par¬ 
ticles. Permutation-invariant momentum-independent 
two-body interactions between particles with arbitrary 
spin-dependence can be decomposed into irreducible 
spherical tensors of ranks 0, 1, and 2 using general re¬ 
lations [34] between irreducible spherical and Cartesian 
tensors. In the present case the Cartesian tensors are 
proportional to the products Sa{j)sa' (f) of the spin com¬ 
ponents of two particles, where a and a' can be either x, 
y, or z. The z-component of the jth particle spin is 

5.(j) = ^(lt(j))(t(j)l-U(j))U(j)l), (1) 

while s„(j) = i[s+(j) -b s_(j)] and Sy{j) = -|[s+(j) - 
s_(j)] are expressed in terms of the spin raising and low¬ 
ering operators 

s+(j) = |t(j))U(j)l, a_(j) = U(j))(t(j)|. 

There are two zero-rank tensors (spherical scalars), 

l>=^y+(r,-r,0 (2) 

i/i' 

and 

=-^(i>. + 14 +_), 

where 

Vzz = “ '^3')Sz{j)sAj') 

j¥=j' 

y+- = ~ 

j¥=j' 


Here 

H±(r) = i(H(r)±H(-r)) 

are the even and odd parts of the two-body potential 
function V{r), 

The three components of the rank 1 spherical tensor 
can be expressed as 

^0^^ = - ^3')s+U)s-if) 

Finally, 

= ^(214, - t>+_) 

^±1^ = =F XI ^^(^3 ~ ^3')Sz{j)s±{j') 

3 ¥^ 3 ' 

^±2 = ^ X ~ I’i')s±(j)s±(/) 

3^3' 

are the five components of the rank 2 spherical tensor 
V^'^K The even or odd rank tensors depend, respectively, 
on the even or odd components H^(r) of the potential. 
Therefore the rank 1 tensor vanish in the case of even 
two-body interaction, while the scalars and rank 2 tensor 
vanish for the odd interaction. 

The interactions conserving the z-projection of the to¬ 
tal many-body spin are expressed in terms of the scalars 
and zero-components of the tensors, 

% ^ X ^(^3 - r,')! t (j))l t (j'))(t (j)l(t (j')l 

i/i' 



% ^ ^ V(r, - ry)| i (j))| i (j'))U (J)IU (j')l 

3 ¥3' 

=+ 5 ^ <■">> 

% ^ X ^(^3 - ry)! t (j))l ; (j'))(t (j)IU (J')l 

3 ¥3' 

I4_ ^ ^ V(ry - ry)| t (j))| i (j'))U (J)l(t (j')l 

3 ¥3' 

= X ^^^^3 - ^rYszU)- 

3 ¥ 3 ' 


( 3 ) 


Here 
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TABLE I. Coefficients in Eq. (10) 


k 


s 



0 

1 1 

2 

0 

3S^-S(S-I-1) 

Sz 

/O S2-S2 

/3[(S-l)2-Sil(S2-S2) 

5(25-1) 

5-1 

\l ^ 5(25-1) 

V 25(5-l)(25-l)(25-3) 

1 


S+2S^+l f 

(S-S^-l)(S-S,) 

/(S-S;,-2)(S-S;j-1)(S2-S2) 


2S{2S-1) 

S-1 V 

2S{2S-1) 

V S(S-l)(2S-l)(2S-3) 

2 

E6(S-S^-1)(S-S^)(S+S;,-M)(S-|-S^-I-2) 

Co 

to 

+ 1)21(S-S;j-2)(S-S^) 

/(S-S^-3)(S-S;,-2)(S-S;j-l)(S-S^) 


25(25-1) 

1 s-iV 

25(25-1) 

V 45(5-l)(25-l)(25-3) 


are zero-components of spherical vectors (see [34]). Equa¬ 
tions (4a) and (4b) lead to the relation 

1 ^ 0 + = ^ (% - ni) • ( 5 ) 


projection, 

(^' 2 '= SsS'SsM<%\y\^us) 



The interactions changing spin of one of the colliding 
particles, 


VTt ^ ^ V{r, - r,OI i (j))l t (/))(t (j)l(t (/)l 
j¥=j' 

^ ^ V{r, - r,,)\ i (j))| i (/))(t Um (/)l 


= -v[\^ - 




( 6 ) 


involve also 


The matrix elements of the spherical vectors and the 
rank 1 tensor follow the same relations (1.23) as the one- 
body interactions. 


{'i’n'S'J^k\'i>„sJ=^S'S.+kXs:k 


\ no 




(9) 


where the factors 


YiS,S',q) 


(_l)S'-S.-fc 


5 y q\ 
-S^-k k ) 

^ [S -S' S'-S J 


^±1 = 

which form a spherical vector together with Vq = + 

Vf (see [34]). 

The interaction 

E__ ^ ^ E(r, - ry)| i (j))| i (/))(t (j)l(t (/)l 

j¥=3' 

= (7) 


changes spins of both colliding particles. Other spin¬ 
changing interactions are obtained by the Hermitian con¬ 
jugation of Eqs. (6) and (7), taking into account that 




( 8 ) 


(S) 

Consider matrix elements between wavefunctions 
with the defined projection Sz of the total spin S. The 
explicit form of 4')^^ , given by Eq. (1.16), is not used 
here. The multi-index n labels different functions with 
the same S and Sz- According to the Wigner-Eckart 
theorem (see [35]), the matrix elements of the spherical 
scalars are diagonal in spins and independent of the spin 


are expressed in terms of the 3j-Wigner symbols. Here 
Ao can be either V^, Vq, or and A±i can be V±i, 
or V^i ■ According to the properties of the 3j-Wigner 
symbols, the matrix elements (9) vanish if [S' — S"| > 1 
(in agreement to the selection rules [19]). The factors 

^S^z’k S' < S are presented in Table 1 in Ref. [1] . 

Hermitian conjugate of Eq. (9), together with Eq. (8) and 
relations (Ho±)t = l/±, (l>o)t = and (Eq^'V = 
give us the matrix elements for S' = S + 1. 

The matrix elements of the components of the rank 
2 spherical tensor E^ can be expressed in terms of the 
matrix elements for the maximal allowed spin projections 
{S' < S) in the same way. 








( 10 ) 


According to the properties of the 3j-Wigner symbols, 
the matrix elements (10) vanish if [S' — S"] >2 (in agree¬ 
ment to the selection rules [19]). The non-vanishing fac- 
tors Xg ’f. , calculated with the 3j-Wigner symbols 

[6, 35], are presented in Table 1. The symmetry proper¬ 
ties of the 3j-Wigner symbols [6, 35] lead to the relation 


y(AS'. 9 ) 

-^Sz-k 


(_l)S-S'-egy(^^AE.9) 
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providing the factors for fc < 0. The matrix 

elements for S' + 1 < S" < S' + 2 are given by Hermitian 
conjugation of Eq. (10), taking into account Eq. (8) and 
the relation ■ 

Thus, each permutation-invariant two-body interac¬ 
tion between particles is expressed in terms of irreducible 
spherical tensors. Their matrix elements for arbitrary to¬ 
tal spin projections are related to ones for the maximal 
allowed spin projections using Wigner-Eckart theorem. 
The next section deals with the later matrix elements. 


II. SUM RULES 

A. Matrix elements for zero-range interactions 

The sums of the matrix elements and sums of their 
squared moduli will be evaluated here for zero-range spin- 
dependent two-body interactions with the even potential 
function 

U(r)=(j(r), (11) 

where r is a H-dimensional vector. This function is gen¬ 
erally used for description of interactions of cold atoms 
in free space, when D = 3, and under tight pancake- 
or cigar-shape confinement, when D = 2 or U = 1, re¬ 
spectively. In the three- and two-dimensional cases the 
i5-function have to be properly renormalized. 

For the even potential function, we have E“''(r) = 
U(r), V~{r) = 0, and, therefore, = Vj,” = 0, 

V+- = VJ^_. Besides, the identity | t (j))(t (j)l + I i 
0))('!' (j)l = 1 leads to the relation 

Vn = l{v-V^-h)- ( 12 ) 

The zero range of interaction allows to relate the ma¬ 
trix elements of Vj-f and _, defined by (4c) and (4d), 

respectively, in the following way. Wavefunctions of in¬ 
distinguishable particles 4* obey to the quantum exclu¬ 
sion principle Vjji'^ = idt , where the sign -|- or — is 
taken for bosons or fermions, respectively. Therefore 

(d/'lU+.ld/) = ± ^(4/'| t (j))l i (/))<5(r, - r,0 

x(; (j)i(t ifmrm. 

The permutation operator Vjj> permutes both spins and 
the coordinates Vj and r^v. However, the J-function sets 
Tj = Tji and the coordinate permutation has no effect. 
Acting to the left, permutation Vjj> of spins is restricted 
by the bra (j, (j)|(t {j')\ of the interaction operator, 
therefore 

(vl/'| 14 _|vl/) = ± ^(4/'| t (j))l i wmrj - r,0 
j¥=3' 

xU (/)|(t (j)l'I')=±(^'|Vtil^)- 


Then (4c) and (4d) for the even potential function, to¬ 
gether with (12), lead to 

(4/'|uJ°)|4/)=-^(4/'|U|vl/) 

(4''|Uo|^')=i(^''|VTT-nil^) 

for bosons. Thus, matrix elements of scalars and zero- 
components of the vector and tensors are expressed in 
term of V, V^i, and Vn- The non-zero components of the 
vector and tensor are obtained from (6) and (7) 

i {V-t - V-i) 

= + ( 14 ) 

= iu__ 

and their Hermitian conjugates. 

For fermions, matrix elements of Uq, and are equal 
to zero, in agreement with the Pauli principle — two par¬ 
ticles with the same spin cannot have equal coordinates. 
This leads to the zero matrix elements of the spherical 
vector and tensors. The matrix elements of the spherical 
scalars are related as 

(vI/'|Uo(°)|4/) = ^(vI;'|U|vl/). (15) 

The sums of the matrix elements of V and sums of their 
squared moduli are derived in Ref. [1]. Other relevant 
interactions are analyzed below. 

B. Matrix elements for non-interacting bosons 

Let us evaluate matrix elements of V^, V^, V-^i, 
and V _between wavefunctions 

= ( 16 ) 

t 

[see Eq. (1.17)] of non-interacting bosons with defined 
projection Sz of the total spin S, where fs is the di¬ 
mension of the respective irreducible representation of 
the symmetric group. The representations and func¬ 
tions whitin the representations, respectively, are labeled 
by the standard Young tableaux r and t of the shape 
A = [iV/2 -I- S', N/2 — S] (see [4, 5]). For bosons, the spa¬ 
tial functions of N non-interacting particles (1.11) are 
represented in the form 
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where V are permutations of N symbols, </Jn(r) are the nal matrices 
spatial orbitals, and the relation for the Young orthogo- 




( 18 ) 


[see (1.7)] is used (see Ref. [1] for other notation). Each 
matrix element, e.g. the one of can be decomposed 
into the spatial and spin parts. 


-r. 0 lCw)( 42 l t (*))l t (*'))(t (*)l(t (*')|sg!)^ 5 .Si. 

t,t' 


Using Eq. (17), the spatial matrix elements can be expressed as 


(19) 


n,Q 


( 20 ) 




The Kronecker 5-symbols appear here due to the orthogonality of the spatial orbitals and the absence of equal 
quantum numbers in each of the sets {n} and {n'}. Due to the 5-symbols, all but two spatial quantum numbers 
remain unchanging. Supposing that the unchanged riin are in the same positions in the sets {n} and {n'}, one can 
see that the Kronecker symbols allow only Q = 7^ or Q = TZVa'. Therefore 




where for the zero-range potential (11) the matrix elements 


( 21 ) 


(n;n2|U|nin2) = (nan'ijUjnins) = / d ripl, {r)(pl, {r)(pn,{r)(p„^{r) 


are invariant over permutations of n( and n^, as well as of ni and n 2 - 

All matrix elements are related by the Wigner-Eckart theorem to ones for the maximal allowed spin projection, 
S'( = S', Sz = S. The spinor matrix elements include projections of the spin wavefunctions (1.14), derived in [36], 




N/2+Sz N 

=c^ss.E^S](^) n it(^j)) n u(^j)), 

•P J = 1 j=N/2+S„ + l 


with the normalization factor 

Css^ = 


{2S + 1){S + Sz)l 


{N/2 + Sz)\{Nl2 - Sy. V (JV/2 + A + 1)(2^)!(5 - ^,)!' 


( 22 ) 


The projection 


Ai Ai N 

(t(*)i(t(*')N^^)=t^5sE^S](^)EV"^upp n it(^j)) n 

V i^v I^j^v i=Ai-ei 

(recall, Ai ,2 = N/2 ± S) can be transformed, using substitution V = QViXi-iPi'Xi, to the form 

Ai Ai-2 N 

(t(*)|(t(*')|Sif^) = C^SsEE^dd](2^'A,-lPcAj5,,2(;,,_l)5.qCA, J] I t (2j)) R U(Qj))- 

Q j=i i=Ai+i 
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The permutations Vi\^-i and Vi'x^ permute symbols in the first row of the Young tableau [0]. Therefore, 

[see Eq. (1.8)] , the summand in the equation above is independent of I and I', 
and the projection can be expressed as 

Ai-2 N 

(t (*)|(t = Ai(Ai - l)C5s5]i?'Jl](Q)d..a(A,-i)5*',QA, n 11 m) n I ^ 

Q i=i i=Ai+i 

The projections involved into matrix elements for other interactions are evaluated in the same way, 

Ai N 

a (*)|(; = A2(A2 - l)CssE^^^I](2)'^^C(Al + l)^^^Q(A,+2) H I t (2j)) H I ^ (^j)) 


J = 1 


i—Ai +3 


Ai-1 


N 


(t(*)ia(*')|Si|’) = AiA2CssE<d](2)^*.SA,5.qQ(A,+i) n lt(2j)) n U(Qj))- 


i=i 


i—Ai +2 


In the spin matrix elements of ■> 


(Sjfjl t (*))! t (*'))(t (*)l(t = [Ai(Ai - l)C55]^E<i(2)<^^C(A.-l)A^QA, 

Q K 

X ^ Sn,QV'V", 
V',V" 

the last Kronecker symbol appears due to orthogonality of the spin states and means that the permutations TZ and Q 
can be different by permutations of particles in the same spin state. They are the permutations V' of the first Ai — 2 
symbols and V" of the last A 2 ones. As the permutations V' and V" do not permute symbols between rows in the 
Young tableau [0], the equality (QP'P") = [see Eq. (1.8)] can be applied, leading to 


(S'fjl t (*))[ t (*'))(t (*)l(t = Ai!A2!Ai(Ai - 

Q 

Let us substitute this equation and Eq. (21) into Eq. (19), perform the summation over t and t', using the relation 


Y,D^^,{V)D^^XQ) = D^Zi'PQ) 


(23) 


[see Eq. (1.6)] and substitute V = Q ^TZ j = TZi, and / = TZi'. Then the Kronecker symbols lead to Vj = Q = 
Ai — 1 and Vj' = Q~^i' = Ai. Equations 

(24) 

(see [5]) and (1.8) lead then to dI’:]JTZVu'Q) = Dl^l^ip-^Q-^V^Q) = h'loi(^”'^AiAi-i) = l?'^|o](iP-i). Then 


using Eq. (18) we get 


Wii'V ^ = VviVi' 

'[ 0 ]' 




) = 2ds'sAi!A2!Ai(Ai - , , (Ai - 1, Ai) 


with 




dA'] 




Matrix elements of other operators are calculated in the same way, 

(Eqiqs'lhllEws) = 2<f5'5Ai!A2!A2(A2 - w+ 1’+ 2) 

(EqiqsE-tlEws) = 2<f5'5-lAl!A2!(Ai - 1)(A2 + 1 )Cs5Cs-15-iEEtE}(^1 - 

= 2 d 5 ' 5 -iAi!A 2 !A 2 (A 2 + 1 )^ 5525 - 15-1 + 1 ) 

= 2<f5'S-2Ai!A2!(A2 + 1)(A2 + 2)CssCs-2S-2V]^'^X{^^{\i - 1, Ai), 

I- 


(25) 

(26) 


(27) 


where A' = [7V/2 + S', N/2 - S"]. 


C. Sums of the matrix elements and their squares 
and products 


Eor the matrix elements which are diagonal in the total 
spin and r, one can calculate their sums and write out 
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them in the form 

X ^(n'n',|y|njny) Y[ 
j<j' 

where Va is any two-body interaction, which does not 
change the spin projection. In each term of the sum 
over j and j', only two spatial quantum numbers can be 
changed. For the operators I^T the factors 

y(^)[VTt] = 2 Ai(Ai -1), = 2A2(A2 -1), 

calculated using the Eqs. (25), (27), (22), (26), (23), 
and (18), are proportional to the numbers Ai(Ai — 1) and 
A 2 (A 2 — 1) of particle pairs with spins f and j,, respec¬ 
tively. For the spherical tensor components, the factors 


are calculated with Eq. (13), 

= \l^s{2S - 1), [vb] = 2S'(iv -1). 

(28b) 

The sum of the matrix elements of the spin-independent 
interactions (1.45a) can be expressed for bosons and the 
zero-range potentials (11) in the form (28a) too with 

y(^) [E] = ^N{N -2) + 2S{S + 1). (28c) 


The sums of squared moduli of the matrix elements 
(25) and (27) and their products are propotional to the 
sums of products of the functions (26), which can be 
expressed as 


E 

r,r' 


n=^3'i 


3\¥^3 i ¥^31 32^32 ¥32 


(29) 


where 




Sg,. (^1, ^2, /^) = E E E D^^[l'iQ)Dllim.,Q3.h,Q3^ (30) 


/ O' g) (g g) '' 

The sums of squared moduli contain different functions S) , namely S] (Ai — l,Ai,Ai — l,Ai) for Ett, 

JIJ 1 J 2 J 2 3i3\32j2 ' / M 




+ + for Vu, (Ai-l,Ai,Ai-l,Ai) for E_t, A:)f^:7;j;)'(Ai, Ai+ 1, Ai, Ai+ 1) for 


,(S-1,S) 


(S-1,S). 


V-i, and A7g^., (Ai — 1, Ai, Ai — 1, Ai) for E_Calculation of the sums of squared moduli of the matrix elements of 

spherical vectors and tensors with (13) and (14) requires also sums over r and r' of the products of matrix elements. 
The latter sums contain — l,Ai,Ai -I- l,Ai 4-2) for products of the matrix elements of E^^ by E^ and 

T’P7E'^(Ai — 1, Ai, Ai, Ai + 1) for VL-t- by VLi. The sums are calculated in Appendix A. 

3 i 3\3232 ' / I -v- 3 i 3\3232 

The sums of products of matrix elements of E [expressed by (1-44) with the zero-range potential function (11)] by 
Eq or Eq are proportional to the sum 


X (n'ln'jEjnj-injj) 5, 

31 ¥ 3 1 ¥31 


(i, i') 

3l3i323^ ' 




32¥32 ¥32 


n'.,, ,n.// 1 
32 ^2 


(31) 


r 


where 


. .,(;,/') 
3l3i3232 ^ ' 


r,r' V 

x4g(^.2.^) (32) 


is calculated in Appendix B. Equation (31) contains 


S^.^l . (Ai - 1, Ai) and (Ai -f 1, Ai -b 2) for Eq 

and Eq, respectively. The sums of squared moduli of 
the matrix elements and their products are expressed in 
different forms if the set of spatial quantum numbers is 


AS) 



changed ({n} {n'}) or conserved ({n} = {n'}). 


D. Changing set of spatial quantum numbers 

If the sets of spatial quantum numbers {n} and {n'} 
are different by two elements,, the product of Kronecker 
symbols in (29) and (31) does not vanish only if ei¬ 
ther ji = j 2 , j[ = j 2 or ji = j[ = j 2 . Since 
12 ) is independent of particular values of 
j and f (see Appendix A), the sum (29) attains the form 


The factors in the sums of squared moduli and prod¬ 
ucts of matrix elements of other operators are expressed 
in a similar form. Explicit expressions for the factors are 
obtained using the normalization factors (22) and sums 

and calculated in Appendices A and B, 

respectively. For example. 


+ 2 


N^{2S+1) 


2S + 3\ 

N{4S^ -f 852 -f S' - 1) -f 2S(2S3 -f 5S2 


S+1 



nihX) 

' r'{n'}r{n}^ J-' r'{n}r{n'}^ 


Equations (13) and (14) lead to the following factors in 
the sums of squared moduli of the matrix elements of 
spherical tensor components and their products 


i<i' 

with 

I 2 ) = V2)+s^;lp(iiXX, k). 

(33) 

Then for any two-body spin-dependent interactions Va 
and Vb, the sums of squared moduli of the matrix ele¬ 
ments and their products can be written out in the form 






(S') 


g(2S-l) / (lV + 2) 
6(2S-h3) V -S'-f 1 


-4S 


=S(^4S + 


iV2_4 

S+1 


Y^^^^\V^^\v]=Asj^S{2S-l) 
r(^’2)[Eo,V'] = 8S(7V- 1) 


(34b) 

(34c) 

(34d) 

(34e) 

(34f) 


^ ’ j<j' 

^ n 

j'¥=3"¥^3 

with S' < S. Each term in the sum above changes two 
of the spatial quantum numbers, conserving other ones. 
Since S and S' are equal to the spin projections, S' is 
unambiguously deterrnined by the operators Va and Vb, 
such that S' = S for y, El)-, and Vn and S' = S + k ior 
vj,^'^ and Efc- The factors Y^^’^'>[Va,Vb] are expressed in 

( O' C'\ 

terms of the sums A 2 ’ ih,li,h,l 2 )- For example, the 

factor y^] takes the form 


y(s. 2 )r+( 2 ) ^( 2)1 ^ {N + 2){N - 2S + 2)(S - 1) 

^ -1’ 2(S+1) 

(34g) 

y('^’2)[y_i,y_i] = {n - 2){n - 2S + 2 ) (34h) 

y(s.2)jy(2)^ = y'^iN -2S + 2)(S - 1) (34i) 

y(s.2)j^_(2)^y^2)j ^ 1(^_25'+2)(1V-2S + 4) (34j) 

The sum of squared moduli of the matrix elements of the 
spin-independent interactions (1.47a) can be expressed in 
the case of bosons with the zero-range potentials (ll) in 
the form (34a) too with 


y(^-")[yTT,l>TT] = 4 [A!A2 !Ai(Ai - l)^^ ' 

X rf’^)(Ai-l,Ai,Ai-l,Ai), 


y(^’2)[y,y] = 6iv(iv-2) + 8S(s + i). (34k) 

The case of a single changed quantum number will be 
considered elsewhere. 


Similarly, the sum (31) leads to the factor 


E. Conserving set of spatial quantum numbers 


y] = 4 Ai!A2!Ai(Ai - 

Is 


(iV-2)! + A7jf4,(Ai-l,Ai) 


If the set of spatial quantum numbers is unchanged, 
{n} = {n'}, the Kronecker symbols in sums (29) and 
(31) are equal to one for any ji, j 2 , j'l, and j^. Then 
these sums contain 
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E E = ^4(E Vn')" + (^3 - 4^4) E - ^^3 + 21 : 4 ) E ’ (35) 

y#ji j2^j^ j'¥^j=A3" j¥^j' 

where the matrix elements of zero-range interactions ( 11 ) are symmetric over permutations of rij and riji, 

= {njnj>\V\njnj>) = {nj>nj\V\njnj>) = {njnf\V\nj'nj) = j (r)v5„^., (r)p (36) 

and Sjij[j 2 j '2 can be either (/i, ^ 2 , ^ 2 ) O'” (^’^^) t®®® (3*^) (3^)] arbitrary superscripts 

and arguments. These functions depend on relations between their subscripts, rather than the subscript specific values 
(see Appendices A and B). Then A 4 = for ji ^ j 2 ^ j'l and ji ^ ^ j(, A 3 = Sjj'jj" + + 

for f y^jy^ j" + f, and A 2 = Sjj'jf + Sjjijij for j ^ j'. The sum (35) can be further transformed to 

E E =iV(iV-l)[(7V-2)(iV-3)A4 + (iV-2)A3 + A2](l^)2 + iV(iV-l)2(A3-4A4)(AiA)2 

h¥=32 

+ N{N - 1)(A2 - A 3 -b 2A4)((A2F)2. 


Here 


N{N - 1) 


is the average value of the matrix elements (36) and 


N 




(A2y)" = 


i=i 

2 


A- 1 


j'y^j 


(38) 


A(A- 1) 


J2(^33'-{V)f 


3<3' 


measure their average deviations [in consistency with 
(1.48)]. Then the sums of squared moduli of the ma¬ 
trix elements and their products can be written out in 
(37) the form 


{yo^'''''\VaMVr +Yl^'°\VaM^lVr +yF'°\VaA]{^2Vy 


(39a) 

with S' < S. The factors here are calculated using results of Appendices A, B, and C. The first factor can be 

represented for all considered 14 and 14 as 


= <5ss'4"^''^[K]l^^'')[14] (39b) 

[see Eq. (28)]. Since S" — S' is unambiguously determined by 14 and 14 (S' = S for A and S' = S -b fc for ' and 

U)) [14,14] = 0 for fc 0 components of the spherical vectors and tensors, when S ^ S'. Other factors are 

expressed as 


^^(S.0)[i>^(2)^i;.j2)] 


rE°^[Vo,Vd] 


(A - 1)(A - 2S)(A -b 2S -b 2)S(2S - 1) / 
3(A-2)(A-3)(2S-b3) V 

(A - 2S)(A -b 2S -b 2)S(2S - 1) A - 1 
6(A-2)(A-3)(2S-b3) V ^ S-b 1 
(A - 1)(A - 2S)(A -b 2S -b 2)S 

sTl 


3444 - 8 S) (39c) 

■ - Ssl (39d) 

(39e) 
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[2{N- 1){N - 2S){N + 2S' + 2)S'(2S' - 1) 
Vs {N-2){S+l) 


,.(S.o)r,>(2) o /2 (A^ - 1){N - 2S){N + 2S + 2)S{2S - 1) 
n [^0 ’^-2^3-(iV-2)(iV-3)- 


,.(s.o)r,>(2) /2 (iV - 25)(iV + 25 + 2)5(25 - 1) 

i'2 y3 (^_2)(7V_3) 

JlV-l)(iV-25)(lV + 25 + 2)5 

n [^0, V] = 2-^-2- 


(39f) 

(39g) 

(39h) 

(39i) 


y^(S.0)j^^2)^^^2)] 

yf-o)[yW,y(2)] 


(iV-l)(iV-25 + 2)(5-l) 


2(1V+ 1)52 - lV(lV + 3) 


{N - 2){N - 3){S + 1) 
{N - 1)(7V - 2S){N + 2S){N - 25 + 2)(5 - 1) 

2{N - 2){N - 3){S + 1) 


N{N - 1){N - 2S + 2) 
y/2N{N - l)(iV - 25 + 2)(5 - 1) 
{N-2) 

{N - l)2(iV - 25 + 2){N -25 + 4) 
{N-2){N -i) 

{N - l)(iV - 25 + 2)(1V -25 + 4) 
2{N -3) 


(39j) 

(39k) 

(391) 

(39m) 

(39n) 

(39o) 


The vanishing factors are omitted above, namely Vh] = 0 when Va or 14 is 14■ The sum of squared moduli 

of the matrix elements of the spin-independent interactions (1.47b) can be expressed for bosons and the zero-range 
potentials (11) in the form (39a) too with 


y/®’0)[l>,l/] 

yf’°)[i>,y] 


(IV- l)(lV-25)(lV + 25 + 2) 
{N - 3) 

{N - 2S){N + 2S + 2) 


3-45 


5 + 1 
N -2 


31V(iV - 4) - 45(5 + 1) + 12 
2{N -2){N -3) ■ 


(39p) 

(39q) 


They are equal to the factors (1.47c). 


III. MULTIPLEX ENERGIES AND 
CORRELATIONS FOR WEAKLY-INTERACTING 
GASES 


Thus, sums of matrix elements and their squared mod¬ 
uli are expressed in terms of universal factors, which are 
independent of the spatial orbitals, and sums of one-body 
matrix elements (or their squared moduli), which are in¬ 
dependent of many-body spins and the spin-dependence 
of the interaction. The sum rules, combined with the 
spin-projection dependencies (9) and (10), provide infor¬ 
mation on each matrix element for any two-body spin- 
dependent interaction between the particles. 


A. Average multiplet energies and energy widths 

Consider a general spin-dependent two-body interac¬ 
tion, conserving the particle spins, 

Vtot = 2ffTr% + 2^11% + 5t4^t4) (40) 

where the potentials are defined by Eqs. (4a), (4b), and 
(4c) with the zero-range potential function (11). Here the 
interaction strengths ga-a' are proportional to the s-wave 
elastic scattering lengthes aaa' for corresponding pairs of 
spin states. For example, in three-dimensional geometry 
(/crcr' = YKh^Qaa' I'm, where TO is the particle’s mass. The 
factors ^ in Eq. (40) appears due to double-counting of 
the interacting pairs in (4a) and (4b). 
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In the case of weak interaction, the average multiplet 
energy can be evaluated in the zero-order of the degen¬ 
erate perturbation theory, in the same way as in the 
case of spin-independent interactions [see derivation of 
Eq. (1.52)] 

Ess. = (41) 


defined by Eq. (37). Substituting the coefficients X and 
Y from Table 1 in [1] , Table I, and Eq. (28) one gets 


EsSz = 9 


-N{N -2) + S{S -f 1) - -S(S +1)^ 
4 3 g 


+ S,iN - l)^ + 

9 9 J 


(V). (42) 


(in the case of spin-dependent interactions the energies 
depend on the total spin projection Sz)- 

For bosons with zero-range interactions, matrix ele¬ 
ments of the potentials are related to ones of the irre¬ 
ducible spherical tensor components by Eqs. (4a), (4b), 
(4c), and (13), 

(vk'lEqlvl/) 

=(4/'| - + iEjvI/). 

These equations allow to expand the interaction Vtot in 
terms of irreducible spherical tensors, which matrix ele¬ 
ments can be related to the ones for the maximal allowed 
spin projections by Eqs. (9) and (10). Then the sum rules 
(28) lead to 


Here the first two terms in the square brackets provide 
the average multiplet energy (1.52) for spin-independent 
interactions. Spin-dependence of the interactions leads 
the third term, which is independent of the total spin 
projection Sz, as well as to the linear and quadratic in 
Sz shifts (the fourth and fifth terms, respectively). The 
corrections are proportional to the ratios 9 ±/ 9 , which 
are determined by the scattering lengths. For example, 
9+/9 ~ —0.001 and g-/g « —0.049 for ®’^Rb. In this 
case, the two states, generally used in experiments, | f) = 
\F = 2,m/ = —1) and | j.) = |F = 1,to/ = 1), have the 
scattering lengths [37] aq « 95.5a_B, ajx « lOO.das, and 
a-i-j, « 98.0aB, where ob is the Bohr radius. One-body 
spin-dependent interactions with external fields lead only 
to linear shifts [see Eq. (1.54)]. 

The root-mean-square multiplet width can be evalu¬ 
ated in the same way as in the case of spin-independent 
interactions [see derivation of Eq. (1.53)] 


Ess. = ^(^gV(^^[V] + g-X^%^’^V(^)[Vo] 

+ y|5+4!o^’'V(^)[Eo(^)])(E), 

where g = (gn + 9u + 9+ = 9n + 9u - ^9u, 

g- = gxt ~ 9 iIj s-nd the average matrix element (V) is 


(AEss.)'^ 


1 

fs 




- E: 


ss.- 


Expanding the interaction Vtot in terms of irreducible 
spherical tensors, expressing their matrix elements in 
terms of the ones for the maximal allowed spin projec¬ 
tions, and applying the sum rules (39) we get 


(AEss.)^ = i E E] + V/^’°)[Vb, Eo] + ^ 


+ 2g.gX^ff^\^^’°\Vo,V] + 2 JE] + 2^ -g+g_X^s%"''^ 


(A.E)2. 


Here the matrix element deviations (A^E)^ are defined by Eq. (38) and the terms proportional to (E)^ are canceled 
due to relation (39b). The first term in the square brackets gives the width for the spin-independent interactions 
(1.53). Corrections due to spin-dependence of the interactions are proportional to the small parameters g±/g. Leading 
terms in the coefficients before each power of Sz can be obtained using explicit expressions for X from Table I in [1] 
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and Table I, 


{AEssy = 






2 5*+ 1 


3 25'- 




9 ^ 


9 


2^6 


5(25-1) 




2^6 


g 52(25-1) 


-^^(S,0) rVV(2) T> 1 ^ o3 

1^0 ^ ^oj o ^z- 


6 


52(25'-!)- 


■\aC'S',0) rT /(2) t /(2)] ■^+ c 4 

;-^2 r 0 5 ^0 Jo *^2 


(A,y)2. 


(43) 


The consideration above was devoted to the case of 
bosons. For fermions, matrix elements of V'\\, and 
vanish, according to the Pauli principle (see Sec. II A). 
Then, due to Eqs. (4c) and (15) matrix elements of Vtot 
are equal to ones of g^iV. Therefore, the muliplet average 
energies and energy widths are independent of 5^ and 
can be calculated using Eqs. (1.52) and (1.53) for spin- 
independent interactions. 


B. Average correlations 


The probabilities of finding two particles with given 
(either equal or different) spins in the same point, the 
two-body local spin-dependent correlations are expecta¬ 
tion values of operators 


= Siri - r2)|CTi(l))|CT2(2))(cri(l)|(cr2(2)|. (44) 


The spin projection aj can be either f or j,. Due to 
permutation symmetry of the total wavefunctions, the 
expectation values of Pcria -2 are proportional to matrix 
elements of the spin-dependent potentials I4 -jCT 2 with the 
potential function (11). The multiplet-averaged corre¬ 
lations can be evaluated in the same way as the aver¬ 
age multiplet energy (41). However, Eq. (41) already 
contains all necessary information, as, according to the 
Hellmann-Eeinman theorem [38, 39], the correlations are 
proportional to derivatives of the average energy over re¬ 
spective coupling constants 


.(AS.) = J_ \o 

2 o - io -2 Jg r{n}Sz'P‘ 




(S) 

•{n}Sz 


1 4- l5crio-2 d - 
N{N - 1) dg„^„^ 


r 


Then the dependence of the multiplet-averaged correla¬ 
tions 




-h 




N{N - 1) 


(P2(0)) 




-h 




N{N - 1) 


1 


N{N - 1) 


(P2(0)) 


on the total many-body spin 5 and its projection 5z is 
factorized to the universal factors, which are indepen¬ 
dent of the spatial Hamiltonian and occupied spatial or¬ 
bitals. These factors are expressed in terms of the coef¬ 
ficients X and Y. The dependence on the spatial state 
is given by the average two-body density [19], (/92(0)), 
which is independent of 5 and Sz and equal to the av¬ 
erage matrix element {V) (37) for the potential function 
(11), (p2(0)) = (V). Similar factorization was proved [19] 
for spin-independent local correlations of particles with 
arbitrary spins. Substitution the coefficients X and Y 
from Table 1 in [1] , Table I, and Eq. (28) leads to ex¬ 
plicit expressions for the universal factors in terms of 5 
and Sz, 


- 2) + 2(iV - l)5z + 25f) 

-1)5^ + 2^.^) 

(45) 

4f =W0^ - 2) + + 1) - 25z^) 

The local spin-independent correlations [19] are 
multiplet-averaged expectation values of (I(ri — r2) and 
can be calculated with the characters from Table 2 in 
[1] for both bosons and fermions (the signs -I- and —, 
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respectively, below) 


pW(0)= (1± 


= 1 ± 


(P2(0)) 

4S'2 + N^ + AS-4N 
2N{N-1) 


(P2(0)). (46) 


They are related to spin-dependent correlations, 
-{S,Sz) , -(S,Sz) , z^-(S,Sz) 1 1 

P 2 (0 ) = PtT + proved 

in the same way as (12), and are independent of Sz- 
In an alternative description, each particle is charac¬ 
terized by its spin projection and coordinate, and the 
total wavefunction is symmetrized for bosons or antisym¬ 
metrized for fermions over permutations of all particles 
[see Eq. (1.19)], 


N 

^{«}W = 

V j=l 

(47) 

where the factor sgn(P) is the permutation parity for 
fermions and sgn(P) = 1 for bosons. Given total spin 
projection Sz, the set {a} contains = N/2 + Sz spins 
t and iVj, = N/2 — Sz spins I- For these wavefunctions 
the correlations are calculated as expectation values of 
the operators (44) and averaged over all distinct choices 
of A^.|- particles with spin up, leading to 


{N^,N^) l+sgn(T’i2) 
NiN-1) 

iN^,N^) _l + sgn(T’i 2 ) 


Pil N{N-1) 

z{Nt,N^) (P 2 ( 0 )) 

-IV(7V-1)^4^4 


iVt(iVt-l)(P2(0)) 

iv^(iv^-i)(P2(0)). 


The transposition parity sgn(Pi 2 ) = 1 leads to the fac¬ 
tor 2 in and for bosons. For fermions, 

sgn(Pi 2 ) = —1 and correlations of particles with the 
same spins vanish, according to the Pauli exclusion prin¬ 
ciple. Using relations between Sz and N^i, one can see 
that the average correlations of bosons with the same 
spins are the same as for wavefunctions with defined col¬ 
lective spins and individual spin projections, = 


However, average correla¬ 
tions of particles with opposite spins are different. 


{S,Sz) _ , (P2(0)) 

+ 7V(]V - 1) 



The same is valid for average local spin-independent cor¬ 
relations, calculated as a sum of spin-dependent corre¬ 
lations. For the defined individual spin projections we 
have for both bosons and fermions (the signs -|- and —, 
respectively, below) 


which depends on the total spin projection, unlike P 2 ^^ (0) 
[see Eq. (46)]. For fermions, the average correlations of 
particles with opposite spins, being equal to a half spin- 
independent ones, are different for the two kinds of states 
too. Thus the spin-independent correlations, as well as 
the correlations of particles with opposite spins, allow to 
determine the kind of the many-body state. 


CONCLUSIONS 

Matrix elements of spin-dependent two-body interac¬ 
tions [Eqs. (4), (6), and (7)] in the basis with collec¬ 
tive spin and spatial wavefunctions (16) can be calculated 
with group-theoretical methods. These matrix elements 
agree to the selection rules [19]. The interactions can be 
decomposed into irreducible spherical tensors, whose ex¬ 
plicit dependencies on the total spin projection [Eqs. (9) 
and (10)] are obtained using the Wigner-Eckart theorem. 
Analytic expressions are derived for sums of these matrix 
elements (28) and their squared moduli [Eqs. (34) and 
(39)] over wavefunctions with the fixed total spin, its pro¬ 
jection, and the set of spatial orbitals. Dependence on 
the many-body states in these sums is given by the 3j 
Wigner symbols and the universal factors 
yd'S'.o) ^ ydS.o) ^ . These factors are independent 

of details of one-body Hamiltonians and are expressed in 
terms of the total spin and number of particles. The sum 
rules can be applied to the evaluation of changes of the 
spin-multiplet average energies (42) and energy widths 
(43) due to weak spin-dependent interactions. Mutiplet- 
averaged two-body spin-dependent correlations (45), cal¬ 
culated with the sum rules, are factorized to universal 
factors, which are independent of the spatial orbitals, and 
the average density, which is independent of many-body 
spins. The difference between these correlations and ones 
for the many-body states with defined individual spin 
projections allows identification of the many-body state 
kind. Other possible applications of the sum rules in¬ 
clude estimates of the spin-multiplet depletion rates due 
to spin-dependent two-body perturbations. 


Appendix A: Calculation of the sums (30) 

The sums (30) and (32) contain the Young orthogonal 
matrix elements (77.), which have been calculated 

by Goddard [40] in the following way. Each permutation 
77 can be represented as 

'^ex 

fe=i 


Ptt 


Pll 


„ (att.wD 


= 1 ± 


7V2 -h 452 - 27V 
27V(7V- 1) 


(P2(0)), 


where V' are permutations of symbols in the first row of 
the Young tableau [0] (Ai first symbols), V" are permu¬ 
tations of symbols in the second row (A 2 last symbols), 
and transpose symbols between the rows as < Ai 
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and > Ai. Then [40] the matrix element is inversely proportional to the binomial coefficient, 


^[o][o](^) = (-1)^ 


= (-ir 




(A2) 

Using relations (23) and (18) and substitution 71 = 
QT’~^, the sum (30) can be represented in the following 
form 


(^1 > ’ ^2, ^ 2 ) = XI ^[0] [0] (^)^[0] [0] (^) X '^'1 ■Ph '^'' 1 TA ,nrn ,nvf^ 


d 0 ][ 0 ]O''/'^[ 0 ][ 0 ]^ 
n V 


where A = [iV/2 + S,N /2 — S'] and A' = [N /2 + S',N /2 — S']. Only the sums with ji ^ j[ and j2 ^ are used 
here. This implies li ^ l[ and I2 ^ l2- The sum remains unchanged on simultaneous permutation of arguments and 

corresponding subscripts, (Zi,/'i, I2, ^2) = ^2, ^2) = ^2)- 

If ji = 32 and j[ = 32 , 




U'l 


(A3) 


n 


since there are {N — 2)! permutations P such that /i = T^ji and = Pj(. This sum is independent of j and 
j'. Due to the invariance mentioned above, fj 12,12) = ^jfjj\h,li,l2,h) and {l{, 1 1,12,12) = 

fj (li,l'i,12,12)- The identity 12,12) = E%^J(l2,l2,li,l'i) can be proved by the substitution 7 Z = 


33 33 
7^-b 

If j'l = J2) but ji ^ J2, we have 


X ^iuV3i^i', ,vj'^h,nvj25i'^,nvj' = X ^i2,ni'^^h,rji^i,vj2^i{,vj' = {N - 3 )ISi>^^tzi'^ X “ ^“1 

“P I 


V 

Then 


^1’ ^ 2 , Q = X ^[0][0](^)'^[0][0](^)(^ - .7^;i(I - K,nw - K,niO 


■R 


__wl'S'b'S)/,/ ,/\ _ 1 /; ;/ 7 ;/\ 

(TV - l)(iV - 2)^" _/v _ 2'^U'U' ^21 


— bj since I2 I2) is independent of j' and ji ^ j2- Here the sum 

Ujf(/,/') = (TV - 1)! X ^[o][o] (^)^[o][0] 


(A4) 


n 


is calculated in Appendix C. The sum 

Sf'^\h,l[,l2, Q = E^f'j}^^(l,,l[,l2,l!2) + E^^'j^^,(l,,l[,l2, 1 ' 2 ) + E^^';p^^(l,,l[,l2, 1 ' 2 ) + (I,, l[, h, /^) 

in Eq. (35) can be expressed as 

1 


Ef’^\h,l[,l2,l'2) = 


(TV - l)(iV _—)^i"'’"\h,l[,l2,1'2) - J^sf'^\h,l[,l2, Q, 


(A5) 


where 


^P'’^^(^1 ,/'i,^ 2,^2) = + E^^f’^\l[,l2) + E^^f’^\h,l'2) + E^^f’^\l[,l'2) (A6) 

and ^2'^ 12,12) is defined by Eq. (33). 
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Finally, if ji ^ j 2 ^ and ji ^ ^ j 2 , the sum is expressed using the following identity, 


V ii¥=i¥=i'i V 

= (N - 4)! ^ i5/2,ra<5;^,ra'(l “ ^ih){^ - Su[){l - 
1,1' 

Then sf 12 , 12 ) [see Eq. (35)] is independent of ji, J2, j[, and and 




n 


N\{N -4)\ 

Is 


<5aA' — 


[ 0 ][ 0 ] 

1 


(TV - l)(TV-2)(TV-3) 


Ef'^\h,l[,h,l'2) + 


1 


(TV-2)(TV-3) 


Sf’^\h,l[,l2,l'2)- (A7) 


[see Eqs. (L5), (33) and (A6)]. 

Therefore, each sum (30) is expressed in terms of (A3) and (A4). 

Let us at first calculate the sums (A3) for S = S'. If h = I 2 = \i — 1, /^ = /^ = Ai, the Kronecker symbols in Eq. 
(A3) select the permutations of the form (Al) with V' which do not affect li and and li ^ i'j, ^ l'^. Therefore there 
are (Ai — 2)! permutations P' , A 2 ! permutations V", and number of distinct choices of the sets of i'j. and i)] are given 
by the binomial coefficients (^^~^) and {^^) , respectively. Then Eq. (A3) can be transformed as follows, 

' ^ex ' 'Tl-ea;/ ' ' 


nex=0 ^ ex / \ ex/ \ 


Al — 2\ A 2 \ Al 

'^ex 


-2 


N\[N -2)\ 
/sAi(Ai-l) 


1 - 


2 A 2 


4 A 2 


(Al — l)(Ai — A 2 + 3) Ai(Ai — l)(Ai — A 2 + 2)(Ai — A 2 + 3)J 


If Ti = ^2 = Al — I, I'l = I 2 = Al, the permutations 




(AS) 


fc=i 


satisfy the Kronecker symbols if V do not affect li and and li ^ i'j. ^ Since Vi^i^ is a transposition of symbols 
in the first row of the Young tableau [0] , Uex = km and 


(Al - I, Al, Al, Al - 1) = (Al - I, Al, Al - I, Al) 

If Ti = ^2 = Al + I, I'l = ^2 = Al + 2, the proper permutations are given by (AI) with V" which do not affect h 
and I'l and li ^ i'l ^ I'l, and with no additional restrictions to V' and i'j. . There are Ai! permutations V' , (A 2 — 2)! 
permutations V", and number of distinct choices of the sets of T]. and i'l are given by the binomial coefficients 
and , respectively. Then 


(Al + 1, Al + 2, Al + I, Al + 2) 


A 2 

(TV-2)! ^ Ai!(A2 

'klex —0 



N\{N-2)\ r 2 

/sA 2(A2 — 1) Al — A 2 + 3 


If /i = ^2 = Al + 1, I'l = ^2 = Al + 2, the Kronecker symbols are satisfied by the permutations (A8) with the same 
restrictions to V" and i'l as in the previous case. Since Pq/j is a transposition of symbols in the second row of the 
Young tableau [0] , Uex = km and 


^jj'jj' (Al + 1) Al + 2, Al + 2, Al + 1) — ^jj’jj' (Al + I, Ai + 2, Ai + I, Ai + 2). 
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If = Ai — 1, I'l = Ai, ^2 = Ai + 1, and l '2 = \i+ 2, the permutations 

n =Pip.p/'p' n (A9) 

fe=i 

satisfy the Kronecker symbol if V do not affect li and and ^ zj, ^ If I 2 ^ i'l ^ for any k, P/pjand Vi'^i'^are 
additional traspositions between the rows of the Young tableau [0], and Uex = km+ ‘^- Otherwise, if one of i'l is equal 
to I 2 , 'Phh'Pi'^h = + 1- By the Same reason Uex = km + 1 if one of zjl is equal to I 2 . If two of i'l 

are equal to I 2 and Z^j '^ex = km- Then 




A 2 

(fV-2)! ^ (Ai-2)!A2! 




A, V' 

km +2/ \km ~ 1/ \km + 1/ 



2N\{N -2)\ 

/sAi(Ai — l)(Ai — A 2 + 2)(Ai — A 2 + 3) 


This derivation is valid for the case of Zi = Ai — 1, Z^ = Ai, Z 2 = Ai + 2, and Z 2 = Ai + 1 as well, giving 

(^1 - 1 , Ai, Ai + 2, Ai + 1) = (Ai - 1, Ai, Ai + 1, Ai + 2). 

Consider now the case of S" = S' — 1. If h = I 2 = Xi — 1, I'l = I 2 = Ai, the Kronecker symbols in Eq. (A3) select 
the permutations of the form (AI) with V' which do not affect Zi and l[ and h I'l- Then 


Y(f-'’^)(Ai-l,Ai,Ai 


A 2 

l,Ai) = (iV-2)! ^ (Ai 

'^ex —0 



Af!(Af-2)!(A?-AiA2 + Ai-2) 

/sAi(Ai-1)2(Ai-A2 + 2) 


If Zi = Z 2 = Ai — I, l[ = Z 2 = Ai, the permutations (A8) satisfy the Kronecker symbols if V' do not affect Zi and l[ 
and Zi i'f, ^ Z(. Since is an additional transposition between the rows of the Young tableau [0] of the shape 
A' = [Ai — 1, A 2 + 1], Uex = km and = fcm + 1- As a result. 


n'oo' 


(Ai — 1, Ai, Ai, Ai 


1 ) 


-(iV-2)! ^ (Ai 

fcm — 0 



-1 


Af!(Af-2)!(Ai + 2) 
/sAi(Ai-1)2(Ai-A2 + 2)- 


If Zi = Z 2 = Ai, Z( = Z 2 = Ai + 1, the Kronecker symbols in Eq. (A3) are satisfied by the permutations of the form 
(AI) with V' which do not affect Zi , V" which do not affect Z(, zj. ^ Zi, and zJZ ^ Z(. Then 


A 2 

'^ex —0 


1)!(A2 



A^!(A^-2)!(Ai-A2 + 1) 

/sAiA2(Ai — A 2 + 2) 


If Zi = Z 2 = Ai, Z( = Z 2 = Ai + I, the permutations (A8) satisfy the Kronecker symbols if V do not affect Zi, V" do 
not affect Z(, zj. ^ Zi, and i'l ^ Z(. Since Vi^i,^ is an additional transposition between the rows of the Young tableau 
[0] of the shape A, Uex = km + 1 and = km- As a result. 


A 2 


+ 1, Ai + 1, Ai) = -(TV - 2)! ^ (Ai - 1)!(A2 - 1)! 


k-rn —0 


Ai — /A 2 — 1 


Ai 

km T 1 


^ "Ai - E 


N\{N -2)\ 


/sAiA2(Ai — A 2 + 2) 
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If h = Xi — 1, l[ = I 2 = Ai, ^2 = Ai + 1, the proper permutations are given by (A 8 ) with V' which do not affect li 
and l[ and li i'/, l[. Vi^i^ is an additional transposition between the rows of the Young tableaux [0] of the shapes 

A and A' and Uex = n'^x = + 1 i unless = Ai + 1, when Uex = n'^x = Then 


us-i.s) 

n'oi' 


(Ai 


1, Ai, Ai + 1, Ai) 


(iV-2)! ^ (Ai -2)!A2! 

fcm=0 




A 2 - 1 \/ Ai A-VAi-IA”' 

+ V \k 

m + 1 / 



m{N-2)1 

/sAi(Ai — l)(Ai — A 2 + 2) 


If li = Xi — 1, l[ = I 2 = Ai, ?2 = Ai + 1, the Kronecker symbols are satisfied by the permutations (A9) with the same 
restrictions to V and as in the previous case. Since PaiAi-iT’aiAi+i = T’aiAi+iT’ai-iAi+i nex = n'^x = + 1 > 

unless zjl = Ai + 1, when Uex = n'^x = km, SiS in the previous case, and 

- 1, Ai, Ai, Ai + 1) = Ai - 1, Ai, Ai + 1). 

The last relevant case is S" = S' — 2. If li = I 2 = Xi — 1, l[ = I 2 = Xi, the Kronecker symbols in Eq. (A3) select the 
permutations of the form (AI) with V which do not affect li and l[ and li ^ i'f. ^ l[. Then 


^{S-2,S) 


(Ai — 1, Ai, Ai — I, Ai) 


(iV-2)! ^ (Ai 

'^ex —0 



N\{N-2)\ 

/sAi(Ai-l)- 


If h = I 2 = Xi — 1, l[ = I 2 = Ai,the Kronecker symbols in 
are satisfied by the permutations (A 8 ) with the same 
restrictions to V and z). as in the previous case. Now 
Pq /2 is a transposition within the same row of the Young 
tableaux [0] of the shapes A or A' = [Ai — 2 ,A 2 + 2]. 
Therefore 


yjiS-2,S) 
'^ 33 ' 33 ' 


(Ai— 1 , Ai, Ai, Ai—I) 


^ ^(S-2.S) 


(Ai, Ai 


1, Ai, Ai —1) 


Eq. (35). 

If ji = 72 , but ji ^ j '2 ^ j'l, 

l^V'^V V 

= (7V-3)! ^ 


Appendix B: Calculation of the sums (32) 

The sum (32) is expressed using the relations (23) and 
(18) as 

V 

Only the sums with ji ^ and j 2 ^ 72 are used here. 
This implies I ^ V. 

If 7 i = 72 and j[ = or ji = j '2 and ji = j' 2 , there 
are (TV — 2)! permutations V which satisfy the Kronecker 
symbols, and the identity (24) leads to 


Similarly, 


.(LI') = S 


(-5) 

313323 


AI, I') = a, I') = (I, I'). 

' ^ > z 313332 z ’ Z 33i332 ' ’ z 


AS) 


li I = Xi — 1, I' = Xi, for I" < Xi — 2, I and I" are in the 
same row of the Young tableau [0] of the shape A, and 
Z7|pj^[P](Pi/") = 1. Eor I" > Xi, I and I" are in different 

rows of this Young tableau, and -D[o]\o](^“") “ ~1/Ai 
[see Eq. (A2)]. Then 


rg>,;(A, - 1,A.) 


(A-3)! 



Similarly, 


33'33 


Ai,n 




{1,1') 


(A-2)!i7W„](P,0- 


33i332 


(Ai + 1, Ai + 2 ) 


(A - 3)! (A 2 - 3). 


Therefore 

- 1 - ^ 1 ) = + 1 , Ai + 2 ) = (A - 2 )!, 

since I and I' are in the same row of the Young tableau 
[0] of the shape A [see Eq. (I. 8 )[. Then 

Yf^(Ai-l,Ai) = rf^(Ai + l,Ai+2) =2(A-2)! (Bl) 


Then 

Yf ^(Ai - 1, Ai) =4(A - 3)! (Ai - 2 - A2/A1) 
E^ (Ai + 1, Ai + 2) =4(A - 3)! (A2 - 3) 

(B2) 

in Eq. (35). 

If 7 i 7 ^ 72 7 ^ j'l and ji ^ j '2 ^ j'l, we have 
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l^V’^V l^V’^V,V"^l" V l^V’^V l^V'^V,V"^l 




Using the same values of the Young matrix elements as in the previous case, one gets 

rf >(A. - 1, A,) =r<3,„;(A. - 1, A,) = (« - 4)’ 
sf\\i + 1, Ai + 2) (Ai + 1, Ai + 2) = (iV - 4)! [(A 2 - 2)(A2 - 5) + Ai(Ai - 1)] 


(Ai — 2)(Ai — 3) + A2(A2 — 1) — 2(Ai — 2) — 

Ai 


(B3) 


Appendix C: Calculation of the sums (A4) 

The sum (A4), expressed as 

ujf= (A- 1 )!^uM 
n 

(Cl) 

is denoted for consistency with the sum UC ’ ' [see Eq. 
(I.A.l)], which is its special case for I = V = Ai. Since 
■^[o][o](^~^) = ^[o][o](^) the sub¬ 

stitution TZ = TZ~^ transposes arguments I and I'. There¬ 
fore 

= (C2) 

Consider at first the case S = S'. If / = /' < Ai, 
the Kronecker symbol in Eq. (Cl) selects the permuta¬ 
tions of the form (AI) with V which do not affect I and 
i'j. ^ 1. Therefore there are (Ai — 1)! permutations V^ 
A2! permutations V", and number of distinct choices of 
the sets of ij. and i)! are given by the binomial coeffi¬ 
cients and (^^), respectively. Then this sum, 

being independent of the particular values of I, is equal 


to Ujj’^\Xi,Xi) = , calculated in Appendix in 

Ref. [1], 


y(S,S) 


iu) 


N\{N -l)\ 
fsX\ 


Ai 


A2 

Ai — A2 -f 2 


If I ^ I', but I < Ai and T < Ai, the Kronecker symbol 
in Eq. (Cl) selects permutations 

f^m 

n = Vw n , (C3) 

fc=i 

with the same restrictions to V and ij.. Since both 
I and V are in the first row of the Young tableau [0], 

in this case. 

If ^ = ?' > Ai, the Kronecker symbol in Eq. (Cl) selects 
the permutations (AI) with V which do not affect I' 
and i'l ^ I'. There are Ai! permutations V' , (A2 — 1)! 
permutations V", and number of distinct choices of the 
sets of zj. and z'^' are given by the binomial coefficients 
and , respectively. Then Eq. (A4) can be 

represented as. 


y,{S,S) 


{1,1) 


A2 

(A-1)! ^ Ai!(A2 

'^ex —0 



(A- 1 )!((A 2 


A 2 — 1 

1)!)’ E 


—0 


(Ai - nex)'. 
(A2 - riex - 1)! 


A!(A- 1)! 
fsX2 



1 

Ai — A 2 -f 2 


(C4) 


If I ^ I', but I > Ai and V > Ai, the Kronecker symbol in Eq. (Cl) selects permutations (C3) with the same 
restrictions to P" and zjf as in the case of / = I' > Ai. Since both I and I' are in the second row of the Young tableau 
[0], in this case. 

If ^ < Ai < I', the permutations (C3) satisfy the Kronecker symbol in Eq. (Cl) if P' do not affect I and i'k ^ I- If 
I'l ^ V for any fc, Piy is an additional transposition between the rows of the Young tableau [0], and rzea; = km + 
Otherwise, if z'^' = V, since PwPiu'^ = Pw^Pw, n^x = km- Then 


uf'')(l,0 = (A-I)! ^ (Ai-I)!A2!(^' ( 

fcm=0 \ m / \ 


I 


Ai 

km + 1 


-2 


-f 


A2 — 1 

km. 1 


-2 


A!(A- I)! 
fsXi (Ai — A 2 -I- 2) 
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The next case is S" = S' — 1. If / = I' = Ai, the sum was calculated in Appendix in Ref. [1], 


^ _ V-(S-I.S) _ N'.{N - 1)! 


If ? = r = Ai — 1, the Kronecker symbol in Eq. (Cl) selects permutations 




(C5) 


fe=i 


if V' are permutations of the Ai — 2 first symbols, V” are permutations of the A 2 last symbols, lo Ai — I, zj, < Ai — 2, 
and > Ai. The numbers of transpositions nex and between rows of the Young tableaux [0] of the shapes A and 
A', respectively, depend on zq. If Zq = Ai, Uexiio) = Kxi^o) = fcm ■ If zo < Ai - 2, Uexiio) = + 1, 

unless zo = z^ for any k. In the last case, Uexiio) = km and, since = 'PXii’f.'PXii'^, ^ea:(*o) = km- Similarly, 

if Zo > Ai, rZg 2 ,(zo) = km , nex{io) = fcm + Ij unless zq = i'l, when nexiio) = km- Thus for 2km + 1 values of zq 
we have nex{io) = ^ea:(*o) = km, for Xi — 2 — km values nex(jo) = km, n'^xi'^o) = + 1, and for A 2 — km values 

nex{io) = km + 1, zZg 2 ,(zo) = km- Then the sum (A4) is expressed as 


A:jf-'’'^)(Ai-I,Ai-l) 


A 2 

(fV-1)! ^ (Ai 


fcm=0 



- (Ai 




-1 


(A 2 



{2km + 1) 



-1 



-1 


/ Ai Nl{N-l)l 

\km + l) 1 A:™ j J /sAi(Ai-I)2- 


If / = Ai, r = Ai — 1, the Kronecker symbol in Eq. (Cl) selects permutations (C3) if V are permutations of the 
Ai — 1 first symbols, V" are permutations of the A 2 last symbols, z'j, < Ai — 1, and zjl > Ai . Now Uex = km , 
n'gx = fcm + 1 unless zj, = Ai — 1 for any k, when Uex = n'^x = km- Then 


Us-i.s) 

^33 


(Ai, Ai 


1 ) 


(iV-1)! ^ (Ai 


fcm=0 



A general relation can be derived for ^ > Ai + 1 and 
arbitrary when the Kronecker symbol in Eq. (Cl) is 
satisfied by 

n = VvnVVin (C6) 

with arbitrary V € Sjy-i, such that PN = N. Then 
Z1 jYp](7^) = nV), since both I and N are in 

the second row of the Young tableau [0] of the shape A 
[see Eq. (1.8)]. If S" < S' and A'^ < Ai, = 

Zl|pj[YP;/ArP) by the same reason. Using Eq. (23), the 
sum can be expressed as 

r,r' 

X E Dl^0]i'P)Di"l]{P)- 

V€Sn-i 

As P are elements of the subgroup Sn-i of permuta¬ 
tions of iV — I first symbols, a reduction to subgroup (see 


/AiA-' [pi-Y pi "A”'- pi-Ypi-Y”' 

\kmJ \km \ km ) \ km ) \km + y 

m{N-iy. 

""/sAi(Ai-l)- 


[4]) can be used, d\^\v) = where the Young 

tableaux f and t, corresponding to the same Young di¬ 
agram, A, are obtained by the removal of the symbol N 
from the tableaux r and t, respectively. = 0 if 

f and i correspond to different Young diagrams due to 
different placement of the symbol TV in r and t .) The 
summation over V can be then performed using the or¬ 
thogonality relation (see [4, 5]) 


E AlPI(Q)i?W(Q) 

QGSn-1 


m 


fx,-xAN-i) 


5tt'5rr'^\y, 


(C7) 

where fs{N — 1) is the representation dimension for N — 
1 particles. The symbol N is placed in the end of the 
second row in the Young tableau [0]. Therefore, A = 
[Ai, A 2 — I], A' = [Ai,A 2 — 1], and A = A' only if Ai = A(, 
or S' = S". As a result. 


Yf’^)(;,/')cx^xv = 0 

whenever ? > Ai -I- 1 and S' < S. Due to Eq. (C2), 
the sum vanish whenever T > Ai -I- 1 for arbitrary I too. 
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This general relation provides the sums appearing in the 
present calculations 


- 1, Ai + 1) = Ai + 1) 

= A:]f-^’^)(Ai +1, Ai) = +1, Ai +1) = 0. 


Another general relation restricts difference be¬ 
tween S and S'. Equations (23), (C6) , and 

the reduction to subgroup lead to (7^) = 

and = 

Then the sum 


(Cl) contains 

^ 45(p)4J(p)cxEv 

VgSn-1 

[see Eq. (C7)]. Now the symbol N can be placed in the 
end of either row of the Young tableaux. Then both A = 
[Ai, A2 — 1] and A = [Ai — 1, A2] are allowed, as well as two 
similar A'. Therefore, |Ai — A'^l < 1, or [S' —S"] < 1. This 
provides the sums appearing in the present calculations 

- 1,Ai - 1) = - l,Ai) 

= A7jf-"''')(Ai,Ai - 1) = A:jf-"’^)(Ai, Ai) = 0. 
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